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The Thermal Pollution Group at the University of Miami has been 
developing three-dimensional mathematical models for predicting the 
hydrothermal behaviors of bodies of water subjected to a heated effluent. 
Generally speaking, these models can be classified in two categories, 
namely, the free surface models which take into account the variation 
in water surface elevations and rigid-lid models which treat the water 
surface as flat. 

To enable the prospective user to use these models as accurate 
predictive tools, particularly to assess the environmental impact in the 
case of cooling lakes, they have to be calibrated and verified at a num- 
ber of sites. The present volume describes the application of the rigid- 
lid model developed by this group to a rather complicated site, namely. 
Lake Keowee in South Carolina. Lake Keowee is rather unique since 
it is used by a Nuclear Power Plant as a cooling lake as well as two 
other hydroelectric stations which use it as lower and upper ponds. 

This is the final verification of these models concluding a series of 
such verifications made possible by funding and technical assistance 
provided by the National Aeronautics and Space Administration (NASA) 
and the Environmental Protection Agency (EPA). 

This model will eventually be made available to all prospective 
users by NASA and EPA. The present volume together with the 
"Three Dimensional Rigid-Lid Model User's Manual" is intended for 
assisting such future users. 


ABSTRACT 


The Rigid Lid was developed by the University of Miami, Thermal 
Pollution Croup, to predict three-dimensional temperature and velocity 
distributions in lakes. This model was verified at various sites (Lake 
Belews, Biscayne Bay, etc.) and the verification at Lake Keowee was 
the last of these series of verification runs. 

The verification at Lake Keowee included the following phases of 
work. 

1. Selecting the domain of interest, grid systems and comparing 
the preliminary results with archival data. 

2. Obtaining actual ground truth and infrared scanner data both for 
summer and winter. 

3. Using the model to predict the measured data for the above periods 
and comparing the predicted results with the actual data. 

The model results have compared well with measured data. Thus, 
the .model can be used as an effective predictive tool for future sites. 
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SECTION 1 


INTRODUCTION 


BACKGROUND 

Understanding the environmental impact of hot water discharge from 
the condenser of a steam power plant is of considerable importance in 
helping to preserve aquatic life and also the overall efficiency of the 
power plant. In recent years, with the improvement of computing facili- 
ties, a number of different types of mathematical models have been de- 
veloped to predict such effects. These models, if properly calibrated, 
are extremely useful in determining the complete environmental impact 
arising due to normal as well as abnormal operating conditions, which 
may be difficult to realise experimentally. This is particularly important 
with respect to estimating these adverse effects during the planning 
stages, thus, helping to minimize them in the final construction. 

The mathematical models solve the basic fluid flow and energy 
transfer equations subject to certain simplifications and assumptions. 

These may be used to determine the long term heat budget for the cool- 
ing lake as well as determine the detailed velocity and temperature dis- 
tributions within it. The accuracy of predictions using these models is 
extremely important and the only way to verify this is to apply ithese models 
to certain known sites and compare the predicted values with actual 
measured quantities. Measurements may be obtained by measuring the 
water velocities and temperatures with anemometers and thermometers 
as well as by infrared remote sensing photuyraphy. Each of these 
techniques has its own advantages and disadvantages. To get a proper 
representative data base, all of these have to be used simultaneously. 

Other factors which affect the water conditions are the meteorological 
conditions, viz. wind speed, solar radiation, air temperature and humidity. 

OBJECTIVES OF PRESENT WORK 

For the past several years the thermal pollution group at the Univer- 
sity of Miami has developed a number of mathematical models for both long 
and short term simulation. The model considered here is a three-dimen- 
sional rigid-lid model designed for short time prediction of detailed three- 
dimensional velocity and temperature profiles in the region of the thermal 
plume. This model has been applied in the past years to several sites, 
namely. Lake Belews (Lee, Sengupta, Mathavan, 1977). The predictions of the 
model at the above sites compared reasonably well with measured data. 

The results of these runs have been published by / the University of Miami in 


earlier volumes in cooperation with NASA. 

The success of this model at the above sites led to the belief that 
this model could be used in future predictions. It was mainly for this 
reason that a site was chosen which was quite complex in nature and 
which would bring out the prediction capabilities and limitations of the 
model. This would also be of great help to future users of the model 
in order to obtain a proper interpretation of the results when used as 
a predictive tool. 

Lake Keowee in South Carolina provides such a site; hence it was 
chosen for the final verification of the model. A description of the site 
is given in the following section. 

DESCRIPTION OF LAKE KEOWEE SITE 

Lake Keowee is located on the north of the state of South Carolina 
about 4G km west of Greenville. It was made from 1968 through 1971 by 
damming the Little River and Keowee River. At present, it constitutes 
Duke Power Company's Keowee Toxaway Complex. 

The lake has two arms connected by a canal (maximum depth 30.5 m). 
There are three power plants on it, namely, the Oconee Nuclear Station, 
Keowee Hydro Station and Jocassee-Pumped Storage Station. The Oconee 
Nuclear Station is a three-unit steam-electric station with an installed 
capacity of generating 2580 MW. The Oconee Nuclear Station draws in 
condenser cooling water from the lower arm of Lake Keowee and discharges 
the heated effluent to the upper arm of the lake. The intake structure 
for the condenser cooling water allows water from 20 to 27 m depth (full 
pond) to pass through. The discharge structure has openings from 9 to 
12 meters deep through which the CCW returns directly to the lower arm 
of Lake Keowee. 

Lake Jocassee is located north of Lake Keowee and is used as a 
reservoir for Jocassee-pumped storage station. The upper arm of Lake 
Keowee also serves as the lower pond for Jocassee-pumped storage sta- 
tion. The Jocassee Station has reversible turbines with a maximum gene- 
rating flow (into Lake Keowee) of about 820 m 3 /sec and a maximum pump- 
ing flow of about 775 m 3 /sec; the net flow into Lake Keowee from Lake 
Jocassee is about 15.5 m 3 /sec. 

Lake Keowee has a full pond elevation of 243. 8 m above mean sea 
level. At full pond it has a volume of approximately 1.18 x 10^ m 3 , an 
area of 74 km 2 , a mean depth of 15.8 m and a shoreline of about 480 km. 
The outflow from Lake Keowee is through Keowee hydro station. The 
flow through Keowee hydro station varies from 1.4 m 3 7sec (leakage) to 
560 m 3 /sec (during peak load operation). The maximum draw down of 
the lake is 7.6 m. 

A map of Lake Keowee is shown in Figure 1. The above data was 
obtained from Duke Power Company 1976. 
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SECTION 2 


CONCLUSIONS 


The main objective of this work was to prepare a mathematical 
model package to assess environmental impact due to a heated effluent. 
Hence, the final test for the model was at Lake Keowee. The reasonably 
good agreement of the predictions with data was shown in the last sec- 
tion. The shape of the plume, as predicted by the model resembles 
closely but does not correspond exactly to the plumes obtained from the 
IR scanner photographs. An important comparison in this case is the 
area covered by each isotherm which gives an indication of the spread 
of the heated effluent. This area is approximately the same in the mea- 
sured and predicted cases for both summer and winter runs. Another way 
the accuracy of the model was determined was by the root mean square 
deviations of the temperature between predicted and measured values 
over the entire domain. These deviations are shown in Table 14. 

The accuracy of the IR scanner isotherms is about 0. 5°C (from the 
sensitivity of the process). The accuracy of the position of the isotherms 
is within ±0. 5 x (grid space) in this case resulting from the lateral dis- 
tortion in the Digicolor map and the process of transfering this on to 
the computational grid for the purpose of comparison. 

Based on the above results, it may be concluded that the predictions 
made by the model are reasonable beyond doubt when applied to cooling 
lakes or a similar site. Hence, the model may be used as a predictive 
tool for obtaining three-dimensional temperature and velocity profiles in 
the vicinity of a thermal effluent discharge, and the results may be used 
in evaluating the performance of existing or future thermal power plants. 


SECTION 3 


RECOMMENDATIONS 


Various numerical models have been developed to study the effects 
of heated discharge and meteorological conditions on bodies of water. 

Most of these models are one or two dimensional. These models have a 
high computational speed but only give horizontally or vertically averaged 
values of temperatures. 

Three-dimensional models, however, have a much finer resolution 
but they consume larger computer time. The three-dimensional rigid- 
lid model can be used to obtain detailed temperature and velocity distri- 
butions in a domain where surface gravity waves are small compared to 
the depth of the domain. This model, as compared to free-surface mo- 
dels, runs faster since surface gravity waves are eliminated by this 
rigid-lid assumption. 

A proper method of using this model would be to run a one-dimen- 
sional model initially to obtain a rough picture of the temperatures and 
then using this model to obtain a better resolution, the 1-D results being 
used as ambient conditions. 

The following improvements have been suggested for the model. 

1. Since all natural flows are turbulent, proper turbulent closures are 
needed to make the model meaningful. At present, the simplest 
possible closures, namely constant eddy viscosities and eddy diffu- 
sivities, have been used. However, better results may be obtained 
by using a higher order closure. 

2. At present, the model uses uniform horizontal grids and stretched 
vertical grids. Nonuniform horizontal grids could be introduced for 
better resolution near the boundaries. 

3. The program has been written to be run as a batch-job on the com- 
puter. It could be made interactive so as to enable the user to run 
it on a terminal. However, this would require some modifications in 
order to reduce the storage space. 
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SECTION 4 


MATHEMATICAL FORMULATION AND MODEL DESCRIPTION 


CHOICE OF MODEL 

Two three-dimensional models have been developed by the thermal 
pollution group at the University of Miami. The first one is a free-sur- 
face model, which takes into account the variation in surface heights of 
a basin. The second one is the rigid-lid model which assumes the verti- 
cal velocities on the water surface as zero. These models have been de- 
scribed in previous publications of this group. The choice of a particular 
model for simulating flows depends on the nature of the site. 

9 

Lake Keowee is a relatively small, closed basin with maximum depth 
of 30 m. The smaii horizontal numerical grid dimensions demand extremely 
small integration time steps to satisfy the Courant-Lewy-Fredrichs condi- 
tions inherent in free-surface formulations. Since the surface waves are 
small compared to the depth, i.e. h/H<<1, a rigid-lid formulation is suit- 
able. This rigid-lid model allows acceptable accuracy with acceptable 
time step size since the C.L.F. criterion is eliminated. This makes the 
rigid-lid model a rather obvious choice. 

The rigid-lid model has been applied previously to Biscayne Bay and 
Lake Belews sites with reasonable accuracy, as mentioned in the previous 
section. This led to the belief that it would be ideally suited for the 
Keowee site. 

DESCRIPTION OF THE MODEL 

(Portions of the following section have been published earlier by 
this group.) 

The rigid-lid model has the following capabilities: 

1. Predict the wind-driven circulation. 

2. Predict the circulation caused by inflows and outflows to the domain. 

3. Predict the thermal effects in the domain. 

4. Combine the aforementioned processes. 

The model solves equations for fluid flow (momentum and continuity) 
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and heat transfer together with the equation of state in a three-dimen- 
sional domain. Since most geostrophic flows are turbulent, the resulting 
Reynolds stress terms in the governing equations are replaced by eddy 
transport coefficients. The fluid is considered to be an incompressible 
Boussinesq fluid and the hydrostatic assumption has been made. 

The rigid-lid assumption imposes a zero vertical velocity condition 
on the surface without affecting the horizontal velocities. This causes 
the surface pressure to be different from atmospheric which under special 
conditions may be translated into surface gravity wave heights if no lid 
is present. This assumption distorts transient time scales but does not 
affect circulation patterns as has been demonstrated by Berdahl (1970), 
Crowley (1969, 1970), Hag and Lick (1973) and Young, Liggett and 
Gallagher (1976). Since surface gravity waves no longer have to be 
reproduced, computer time is saved considerably and this assumption is 
quite adequate for cooling lake studies where surface waves are not 
large and do not change rapidly. 

To convert regions having irregular bottom topography to constant 
depth regions for computation stability a vertical stretching, suggested 
by Freeman et al. (1972), has been used. 

GOVERNING EQUATIONS 

The vertical stretching used in the rigid-lid model, originally 
adopted by Sengupta and Lick (1974), is of the form 


Y 


Z 

h(x,y) 


This enables the same number of grid points to be used all over 
the domain without using variable grid sizes. The resulting nondimen- 
sional governing equations in the new (a, 0,y) coordinate system instead 
of (x,y,z) are 


Continuity Equation: 
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Momentum Equation: 
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Hydrostatic Equation: 
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Energy Equation: 
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Equation of State: 

fS = 1.029431 - 0.000020 T - 0.0000048 T 2 
where £ is in gm/cc and T is in °C 
and 


where 

and 

where 
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w(2=0) = 0 (rigid lid) 
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Quantities with subscript 'ref 1 are reference quantities; H and L are 
vertical and horizontal length scales. The variables with wavy lines on 
top a ris dimensional quantities. 

If A h = A ref and = B ref then A* = 1 = B* 





-~fT' 


P 


r 




Uref 1 - 


B 


= -aH 


ref 


u 


U 


ref 2 


If Prandtl number is equal to 1, then A ref = B rg ^. 

A^ and A v are the eddy viscosities in horizontal and vertical directions. 

B^ and B v are the eddy diffusivities in horizontal and vertical directions. 

To obtain a predictive equation for pressure, the horizontal momentum 
equations are integrated from Z=0 to Z=h, where h is the nondimensional 
depth h/H. The integrated equations are then differentiated with respect 
to a and 8 summed. 


The Poisson equation for surface pressure becomes: 



3 2 P 

TP? 


- J. A 

h3a l X 


+ A + C 


1 


- V 


+ 1 — r-A 
h 9 B' y 1 


- A 


1 rill 3h 3P s 

K 1 3o 3a 3 3 3 8 


+ C - Y 

y 

n 9 (ft) 

1 3 1 


(Z=0* 


The last term is the Hirt and Harlow (1964) correction term which accounts 
for nonzero vertical velocities at the rigid lid. The variables (B , B , 
and A , A etc.) are given below: x Y 

X 1 x 2 

A x, = &' { l^ huu) + ll> huv > + h f^ nu > Wf 
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c x = K- i'k < h f 0 ) + h ih Tp f f> Wy 

e 

X p " E uo' 1h{ fT ^ + 4^* - t&pWy 

A = ^ ^fa^ huV ^ + foJ huV + h " 3 ^ T ^ y 

y t o 3 a da oy 

A = /' udy 

y 2 B o 


S = r- / + fr*i§» + £ i id Y 


3U, 
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r 3U, 


Y p ■ E UO /1h{ to /YpdY + h !so /YpdY - Y t PldY 


B - = E ufa ^ TpdY + E u h f? /YpdY " E ” Y ^ p 


3h 

•u’3? 


The set of equations (24-30) together with appropriate boundary 
conditions constitute the mathematical model. The model has been de- 
scribed in earlier publications by this group (Lee and Sengupta, 1976). 

INITIAL AND BOUNDARY CONDITIONS 


The nature of the equations requires initial and boundary conditions 
to be specified. As the initial condition, the velocities, temperatures 
and densities are specified throughout the domain. Boundary conditions 
are specified at the air-water interface, geographical boundaries of the 
domain, the bottom of the basin and the efflux points. At the air-water 
interface a wind stress and a heat transfer coefficient are specified. 

The conditions on the lateral walls are no-slip and no- normal velocity 
for the momentum equations. These walls are assumed to be adiabatic. 

At the floor of the basin, the conditions of no-slip and no-normal velo- 
cities are applicable. The energy equation has a heat flux boundary 
condition. 

At points of efflux open-boundary conditions are specified. If the 
flows are known at the points of efflux, the flow velocities are specified. 
Usually, the temperature is known only at the point where the heated 
discharge (from the power plant) enters the domain. At points where 
the temperatures are not explicitly known, open-boundary condition (i.e. 
zero gradient condition) is specified. The same holds true if the velo- 
cities are not explicitly known (e.g. the connecting canal between the 
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two arms of Lake Koewee) . 

Hence, the boundary conditions are: 
At the surface 


y = 0 

ft = 0 (Rigid Lid) 


3U _ . 
3y * 1 


hH , 

IT )t zx 
ref v 


hHK 

3T _ , s 
3y ~ ' B 


) (T e - T s ) 


3 v r hH . , 

U ref Av Tzy 


where t 2x and are wind stresses in the x and y-directior 

Tg is the equilibrium temperature. 

T is the water surface temperature, 
s 

K g is the surface heat transfer coefficient. 

t . and x are calculated from the wind velocity usin< 
zx zy curve (B. W. Wilson, 1960). 

K and T are calculated as follows: 
s e 

K = 4.5 + 0.05 T + BF(w) + 0. 47(w) (wm" 2 °C _ 
s s 

where T is in °C. 
s 

F(w) = 9. 2 + 0. 46 W 2 (wm 2 mmHg 1 ). 

B = 0.35 + 0.015T + 0.0012 T 2 (mmHg/°C). 

m m 

T m = (T s +T d )/2 <° C) - 

= Dewpoint temperature (°C). 

E d s s 

H = Surface solar radiation (w/m 2 ). 
s 


respectively. 


Wilson 


) 
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At the bottom of the basin 

y = 1 

n = o 

u = 0 

v = 0 

3T =0 
3Y 

On lateral walls 

u = 0 
v = 0 

£ 2=0 

3T __ 3T y 3h 9T _ , 

3 X ~ 3 a h 3 a 3 y 

3T_ _ 3T_ y 3h 3T 
3y ‘ 88 ' h 3MY " 


SPATIAL DIFFERENCE SCHEMES 

The numerical solution of the momentum and energy equations are 
explicit. The values of velocities and temperatures at a future time are 
determined completely using the values at the present and previous 
time steps. Finite difference forward time and central space difference 
schemes are used. Diffusion terms are written using a Dufort-Frankel 
finite difference scheme to relax the diffusive stability criterion. The 
convective stability criterion, however, is not affected. 

A horizontal staggering is used in the computational grids. Hori- 
zontal velocities and temperatures are calculated at the main- grid points 
while vertical velocities and pressures are calculated at half-grid points 

The predictive Poisson equations for calculating rigid-lid pressures 
is finite differenced using a five-point scheme. This is solved by 
successive over relaxation (Liebmann Method). Terms on the right hand 
side of the pressure equation are obtained by integrating terms in the 
horizontal momentum equations over the depth using the trapezoidal rule. 

At the boundaries, single-sided schemes are used since the boundary 
points do not have two adjacent points. A curve is fitted through the 


two most adjacent points towards the interior of the domain. The values 
of the variables, where they are not defined, are obtained by averaging 
the values at four points around the point where the values are defined. 

STABILITY 

It is not possible to make a strict stability analysis of the system 
of equations under consideration. It is customary, however, to take 
advantage of the stability analysis for the one-dimensional Burger's 
equation since this contains an unsteady terrr, a convective term and a 
diffusion term. 

In the present case these criteria can be written as follows for the 
choice of the time step A? 

Convective . - Ax 

At < ~\T 

where U is the maximum horizontal velocity in the domain. 

Diffussive At < 2 


MARKER MATRICES 


Since natural bodies of watjir have irregular boundaries, the physical 
boundaries have to be approximated in a rectangular coordinate system 
using marker matrices (Fortran variable MAR for the full grids and MRH 
for the half grids). The convention used is as follows: 


MAR = 0, point outside the region of interest (i.e. , on dry land) 


In approximating such a boundary using a rectangular grid system 
only a portion of the resulting grid falls within the water. To prevent 
computations to be carried out on dry land markers have to be used to 
distinguish between points lying within and outside the region of in- 
terests. Such markers have to be used both for the main and half-grid 
points. Fortran symbols used are MAR for the main grid and MRH for 
the half-grid system. The invention used is as follows: 


MAR = 0, 
MAR = 1, 
MAR = 2, 
MAR = 3, 
MAR = 4, 


point 

point 

point 

point 

point 


outside the region of interest 
on far y-boundary. 
on near y-boundary. 
on near x-boundary. 
on far x-boundary. 


(i.e., on dry land). 
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MAR 

MAR 

MAR 

MAR 

MAR 

MAR 

MAR 

MRH 

MRH 

MRH 

MRH 

MRH 

MRH 

MRH 


= 5, outside corner, on near x-boundary and far y-boundary. 
= 6, inside corner on far x-boundary and far y-boundary. 

= 7, outside corner on near x-boundary and near y-boundary. 
= 8, inside corner on near x-boundary and near y-boundary. 

= 9, outside corner on far x-boundary and near y-boundary. 

= 10, outside corner on far x-boundary and far y-boundary. 

= 11, interior points (within region of interest). 

Similarly for the half-grid points. 

= 1, corner on far x-boundary and far y-boundary. 

= 2, points on near y-boundary. 

= 3, points on near x-boundary. 

= 4, corner at near x- and near y-boundaries. 

= 6, far corner on x-axis. 

= 7, corner at far x- and y-boundaries. 

= 9, interior grid points. 




SECTION S 


APPLICATION TO LAKE KEOWEE 


INTRODUCTION 

The three-dimensional rigid-lid model described in Chapter II has 
been applied to Lake Keowee, South Carolina, to predict the wind-driven 
circulation, circulation caused by inflows and outflows, and the thermal 
dynamics of a region of this lake. This region includes the buoyant 
plume, the Keowee hydro dam, the connecting canal and an open-boun- 
dary where the effects of inflow and outflow to Jocassee-pumped storage 
station are felt. 

Lake Keowee is a warm monomictic lake having one circulation period 
during the year, beginning in fall, and is thermally stratified during the 
summer. During the circulation period, the water mass is vertically 
mixed. This type of lake can also be termed holomictic (Hutchinson, 

1957). Lake Keowee could also be classified as a subtropical lake as 
it exhibits thermal stratification, a period of total circulation preceded 
by the fall overturn, and surface temperatures above 4°C. The highest 
surface temperatures typically occurred in July and August and ranged 
from approximately 27 to 29°C. During the periods from November through 
February the lake was isothermal or nearly so. It can be concluded from 
the above that a series of meteorological events primarily govern Lake 
Keowee's thermal characteristic. However, as will be shown below, the 
near field region which makes up the region of interest in this present 
study, is very strongly affected by the buoyant plume, inflow and out- 
flow and the Jocassee-pumped storage station. 

CHOICE OF DOMAIN AND GRID SYSTEM 

The region of interest has already been mentioned in the last sec- 
tion. Here, an attempt will be made to justify this choice and then a 
description of the grid system will be given. 

The study and prediction of the thermal impact on the hydro-dyna- 
mics of Lake Keowee has been done in two main sections: 

T. Using a 1-D continuous temperature profile prediction model over 
the period 1971-1979, the area of interest includes the two main 
arms of Lake Keowee and the connecting canal. The application of 
this model has been published by the thermal pollution group of the 
University of Miami (Sengupta et al 1980). 


14 


2. Using the 3-D rigid-lid model, which is described in this report. 

The division has been necessary since the 1-D model predicts 
temperature profiles continuously over a period of time (several years 
if required) ; it can be used to predict the initial conditions required to 
run the 3-D model. The 3-D model, on the other hand, has better 
resolution and is therefore used for predicting the near-fieid thermal 
impact on the lake. This near-field area of interest is shown in Figure 
5. In selecting this area some preliminary runs were made. These are 
described later. 

The region of interest is divided into a square grid as shown in 
Figure 5. The 1 and J axes (corresponding to u and v-veiocity axes) 
are numbered as shown in the above figure. For the top 'open boundary', 
'J' increases from 7 to 20. The left side boundary shown '!' increasing 
from 1 to 17; the right side boundary shows M' increasing from 1 to 16, 
and finally, the bottom boundary shows 'J' increasing from 1 to 18. The 
orientation of these boundaries with the north-south direction is also 
shown in this figure. The squares are (152.4 x 152. 4) square meters. 
Arrows indicate directions of flow into or out of domain. 

SUMMARY OF DATA 

For the purposes of calibrating and verifying the model, an archival 
data base was established and two remote sensing and field data collection 
missions (summer and winter) were undertaken. 

Archivai Data Base 

The archival data base are summarized in Figure 6 and Tables 1 
through 4. These figures were taken from Duke Power Company (1976). 
Figure 6 shows the measured surface isotherms for September 10, 1975. 
This figure constitutes the data base on which the comparisons for the 
preliminary and archival runs are based. The runs are described below. 
Table 1 shows the monthly gross thermal capacity factors for Oconee 
Nuclear Station (1973-1 977). Table 2 is a summary of Oconee Nuclear 
Power Plant for September 10, 1 975. Table 3 shows the input data 
used for the preliminary and archival runs. The data include the flow- 
rate, inlet temperature, discharge temperature and velocity, ambient 
temperature, depth at discharge, discharge width, air temperature, 
wind speeds and the vertical and horizontal eddy diffusivities. Table 4 
is a summary of the volume and area data for Lake Keowee. 

August 1978 Data 

To obtain an understanding of seasonal behavior of Lake Keowee, 
data was gathered both in summer and winter. This data was used to 
verify temperatures predicted by the mode!. Two extreme weather 
conditions were selected with the idea that if the mode! predictions were 
accurate for these conditions it woiild be so for other intermediate con- 


ditions (viz., spring or fall). 

In spite of the fact that the model was used for predictions in the 
vicinity of the thermal plume, temperature and velocity readings were 
taken for the entire lake. This was done in order to ensure complete- 
ness of the data and also to check the choice of the domain of interest 
within the lake. 

Summer data missions were carried out during the period of August 
24 and 25, 1 978. The data collection team consisted of representatives 
from NASA, EPA, University of Miami and Duke Power Company. Three 
boats provided by Duke Power Company were used in the collection of 
ground truth data. The data collection stations are shown in Figure 7. 
This included measurements of water temperatures at various depths 
(from the surface to about 30 m) using YSI type thermistors and velocity 
measurements using an Endeco Type 110 current meter. The thermistors 
were calibrated before each set of readings using a mercury thermometer. 
At each measuring station the boat was anchored and the thermistors 
and velocity meters were lowered by cables. The cables were marked 
thus indicating the depth below the water surface where readings were 
taken. The current meter indicated both the magnitude and direction 
of the horizontal component of the water velocity. The main problem 
encountered here was that most of the velocities were very small and 
close to the threshold of the instrument. Hence, drogues were used to 
determine average surface velocities. 

While the ground truth data was being collected, overflights were 
made by NASA aircraft to obtain synoptic IR scanner data. The aircraft, 
a twin engined (Beechcraft) (NASA-6), was specially equipped for this 
purpose. Flights were made at altitudes of 1000, 2000 and 3000 ft. A 
diagram of the flight plan is shown in Figure 8. The scanner photo- 
graphs were taken with window settings between 78 and 84°F. To cor- 
rect the IR scanner photographs with the ground truth data for each 
set of flights, the water surface temperature was measured with the 
aircraft overhead for at least one station. This was to correct the scan- 
ner readings due to errors caused by water vapor attenuation. 

The scanner used was a daedalus series DS-1250 and remote sensing 
of 8-14 vim radiation is achieved by a Hg: Cd: detector. The detector 
had a 0.015 inch square sensitive area, which was optimum for the reso- 
lution and temperature sensitivity required. This detector was mounted 
in an end-looking, metal cored dewar which had sufficient liquid nitrogen 
coolant capacity for approximately six hours of operation before refilling. 
This system projected through the floor of the airplane, NASA 6, and 
had a scan angle of 120° centered above the vertical. A horizontally- 
mounted telescope with its axis along the direction of flight of the air- 
plane was contained within the scanner. A mirror rotating at 3600 RPM 
and mounted at 45° to the telescope directed heat radiation from the 
ground into the system. A one-third revolution of the mirror covered a 
complete step perpendicular to the scanner axis. Optical resolution 


obtained by this method was about 1.7 milliradians, so the ground areas 
detected become a function of flight altitude; the data accuracy is 0.5°C. 

The video signal from the infrared detector was amplified and re- 
corded on magnetic tape in the aircraft. A method developed by Daedalus 
called Digicolor was used to convert this stored information directly into 
color coded strip imagery. This process limits the number of output 
colors to light (from white 'hottest* to black 'coldest') for any input con- 
dition. The six colors between white and black indicate the six calibrated 
levels of the set of interest. The scanner's thermal reference sources 
were present in flight to 66°F and 84°F respectively, and the settings 
were recorded on the same track with the detector video to insure accu- 
rate voltage relationships irrespective of all amplifier gain adjustments. 

The color bands in the final Digicolor map indicate zones of con- 
stant temperature within the accuracy of resolution of the scanner; hence, 
the line formed by the junction of two adjacent color bands indicate an 
isotherm. The actual temperature of the isotherm is obtained by adding 
a ground truth correction term to the temperature indicated in the map. 

Altogether, three IR flights were made during the summer data 
collection mission. Out of these, two runs were made on August 24. 

The first run was from 0853 hrs EST, On August 25, a single run was 
made from 0908 hrs to 0953 hrs EST. 

The Digicolor maps were transfered on to enlarged maps of the region 
of interest so as to obtain a map of surface isotherms which could be com- 
pared directly with the values predicted by the model. Out of the three 
runs made, only two, namely August 24 morning and August 25 morniing, 
were used for comparisons since the resolution of the colors in the re- 
maining Digicolor map was very poor. 

The values of Oconee Nuclear Stations flows and temperatures every 
hour are obtained from continuous water quality monitoring stations of 
Duke Power Company. Flow through Jocassee-pumped storage station 
and Keowee Hydro Station as well as meteorological conditions (viz. , air 
temperature, wind speed and direction, humidity and incident solar radi- 
ation) were also obtained from continuous monitoring stations. The flows 
through Jocasse, Keowee and Oconee as well as the discharge temperatures 
are shown in Tables 5 and 6. The meteorological data (obtained hourly) 
are shown in Tables 7 and 8. Figures 9 and 10 show the hourly variation 
of Keowee Hydro Station flowrate on these two days. Figures 11 and 12 
show the same for Jocassee-pumped storage station. 

February 1979 Data 

The winter data mission was carried out during February 27, 28, 

1979. The same equipment were used in this mission. Three boats 
were used to measure water temperatures and velocities up to a depth of 
30 m. Two of the boats were equipped with thermistors only and were 
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used to measure temperatures in both branches of Lake Keowee. The 
third boat carried the current meter as weli as a thermistor. The stations 
where readings were taken are indicated in Figure 13. Station #13 was 
used for ground truth correction of the IR scanner data. This point was 
chosen since it is virtually unaffected by the discharge, and the tempera- 
tures here remained fairly constant. 

Ground truth measurements were taken from 0950 hrs to 1353 hrs 
EST and 1600 hrs to 1833 hrs EST on February 27. Due to technical 
problems with the aircraft (NASA- 6) the IR flights were delayed and 
were only from 1549 hrs to 1711 hrs EST on this date. The flight plan 
was identical to that used in the August data, mission. Black body set- 
tings of 38°F and 74°F were used and the flight altitudes were 2000, 

3000 and 1000 feet. 

On Feburary 28, 1979 ground truth measurements using three boats 
were taken from 0851 hrs to 1201 hrs EST and 1420 hrs to 1830 hrs EST. 
IR flights were run from 0850 hrs to 1002 hrs at altitudes of 2000 and 
3000 ft. The black body settings used were established by NASA. IR 
isotherms were constructed for the domain of interest at the University 
of Miami. These maps were used for verifying the results predicted by 
the computer. 

For obtaining meteorological data and flows through Oconee Nuclear 
Station, Jocas see- pumped storage station and Keowee Hydro Station, 
continuous monitoring stations of Duke Power Company, were used. 

Since the velocities in the lake were found to be extremely small, drogues 
were also used. The flows through the three power stations are shown 
in Tables 9 and 10. The data obtained was hourly and the variation of 
the flows through Jocassee and Keowee are shown in Figures 14 and 15 
and Figures 16 and 17 respectively. The meteorological data collected on 
February 27 and 28 are shown in Tables 11 and 12 respectively. 

The ground truth data collected during the summer and winter 
missions are presented by Duke Power Company (1978 and 1979). 

CALCULATION OF INPUT 

The 3-D rigid-lid model as described in Section 4 solves the three- 
dimensional momentum, continuity and energy equations. As shown be- 
fore, these equations are a set of nonlinear, coupled partial differential 
equations and require the following for complete solution. 

1. Initial values of the velocities and temperatures must be specified 
at ail points within the domain. 

2. Boundary conditions for the above variables must be specified at 
all boundaries. 

The choice of the domain of interest as well as the grid system has 
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been described in earlier sections. At ail solid boundaries the velocities 
and temperature gradients are specified as zero. At the surface the 
rigid-lid constraint of zero vertical velocities is assumed. The initial 
conditions assumed for starting the runs are zero velocities and con- 
stant temperatures everywhere within the domain. For subsequent runs 
the results of the previous run is used as the initial condition. Hence, 
the quantities yet to be specified are: 

1. Temperatures and flow velocities at the Oconee Nuclear Station dis- 
charge. 

t 

2. Temperatures and velocities at the Keowee Hydro Station. 

3. The same at the Jocassee boundary and at the canal. (Items 1 and 2). 

4. Surface horizontal velocities and temperatures. 

The above quantities are termed as inputs to the model and the 
procedure used for obtaining them are discussed briefly below for the 
sake of completeness. For further details regarding the actual running 
of the programs, the reader is. advised to refer to the 3-D riqid-lid 
User's Manual. (Sengupta et al, 1980) prepared for this purpose. 

The calculations shown below were for February 27, 1979 simulations. 

Reference Quantities Used 

Reference length = L = maximum length of the domain = 2895.6 m. 

4/3 

Reference horizontal eddy viscosity A ref = 0.002 L s= 38311.48 cm 2 /sec. 

(Note: The constant '0.002' changes with different sites. This particu- 

lar value was used in running the model at Lake Belews and at Biscayne 
Bay. In this case the best value of the constant was found to be '.003' 
which yielded a value of A rgf = 60,000 cm 2 /sec. 

Reference depth = H = max depth considered = 16 m 

Reference vertical eddy viscosity A = 0.002 x H = 37.43 cm 2 /sec. 

Reference velocity U ref = 30 cm /sec. 

Reference time T rg ^ = L/V ref = 9652 sec. 

Oconee Nuclear Station Discharge Velocity 

The discharge is considered to take place through a point at a depth 
of 12 m (k=3). The discharge velocity is calculated as follows. 
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■1 





152.4 m 


The total discharge = {100 x V(— — ) x 152.4 x 12) = Q 

m sec 

(where Q = average discharge in m 3 /sec.) 

. * .V = 7. 4207 cm /sec 

The average value of 'Q* over 24 hrs is taken since the variation 
is negligible. 

V V 

Nondimensionai discharge velocity ~ ^ = -^q = 0. 24740. 

ref 

Keowee Hydro Discharge Velocity 

The outflow through the Keowee hydro station is through a channel 
(152.4 m) x (12 m). 

The volume flowrate Q = (152. 4 x 12 x V) m 3 /sec. 

(where V = discharge velocity (m/sec.) 

V = [0/(152.4 x 12)] m/sec = 152 . 4 .° 2x10 (, cm /sec 

Q is specified as a function of time with the help of polynomials and 
other functions. The curve is approximated and specified in subrou- 
tine INLET1; (refer to the user's manual (Sengupta et al, 1980). 

Keowee hydro flow approximation: refer to the user's manual (Sengupta et al , 1 980 ) 

(February 27 Data) 

SX = Conversion factor for converting discharge (incfs) to nondimensionai 
velocity. 

SV1 = Nondimensionai velocity. - - - - - - - • 

The velocities are approximated as follows: 

SV1 = 0.048 0 < TSDT < 6 
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SV1 = SX * (((17. 54 - 0.048)72.) * (TSDT - 6.0) + 0.48) 

6 < TSDT < 8 

SV1 = SC * (((.048 - 17.54)74.) * (TSDT - 8.0) + 17.54) 

8 < TSDT < 12 

SV1 = SX * (0.048) 12 < TSDT < 24 

Jocassee Flow Velocity 


The entire flow to or from the Jocassee- pumped storage station is 
assumed to take place through the entire upper boundary. The flow 
through this area is shown in the following figure. 

T 

1 6 m 

Jl. 


The flow is assumed to be uniform over this area (i.e., equal flow 
velocities at all internal grid points within this area) and is assumed to 
take place simultaneously with the flow through Jocassee -pumped storage 
station. 


V 




«13 x 152. 4 m 

> 


V = Q /[ (16 x 13 x 152. 4) x 100] cm /sec 

(Q = flow through Jocassee- pumped storage station (m 3 /sec).) 

Q is positive when Jocassee is generating (i.e. # the flow is into the 
region of interest) and negative when pumping (i.e. , flow out of the 
region of interest) . 

Jocassee flow approximation: 

(February 27 Data) 

TSDT = Time from start of run (hrs) . 

SV = Velocity of flow through Jocassee boundary (nondimensional) . 

SF = Conversion factor to convert flowrate (cfs) to nondimensional 
velocity 

= 0.00322579 

The velocity is approximated as follows. 


SV a SF * (-14.395 - (18.75 - 14.395) * (TSDT) 

1 < TSDT < 0 

SV = SF * (-18.754) 1 < TSDT < 5 

SV = SF * (((16.823 + T8. 75,4)73.).,* (TSDT t 5.Q) + 18.7.54) 

5 < TSDT < 8 

SV = SF * (((16.823 - 0.1) /3.) * (TSDT - 8.0) - 16.823) 

8 < TSDT < 11 

SV = SF * 0.1 11 < TSDT < 23 

SV = -SF * ((4.5 + 0.1) * (TSDT - 0.1) 

23 < TSDT < 24 

A similar procedure was followed for simulations of the other days. 
Condition at Open Boundaries 

Open boundaries are those where the values of temperatures and/or 
velocities cannot be specifically obtained but continuity of flow has to be 
maintained. One such boundary is at the mouth of the canal connecting 
the two arms of Lake Keowee. The condition — = 0 is sepcified here 
both for velocities and temperatures. At the Jocassee and Keowee boun- 
daries, the same kind of zero gradient conditions is specified for tempera 
tures only. The calculation of parameters (T-, TAU, etc.) for the sur- 
face boundary conditions are shown in Section 2. 
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SECTION 6 


RESULTS AND DISCUSSIONS 


The results of simulation using the 3-D rigid model at Lake Belews 
and Biscayne Bay are shown in previous publications by the University 
of Miami (Lee, Sengupta and Mathavan, 1977, and Lee and Sengupta, 
1976). In the above two cases the predictions made by the model agreed 
closely with IR scanner and ground truth data. The following section 
discusses the runs made for Lake Keowee. 

Preliminary Runs and Results 

These runs were primarily made for the selection of the boundary 
conditions, initial conditions and to test the behavior of the model when 
incomplete and/or arbitrary data are used. These runs are summarized 
in Table 13. 

1. Run Number L001 

(The Number L001 is a sequence number and must only be interpreted 
as such.) 

This run was essentially made for debugging the computer program 
modifications. The features include a discharge velocity of 5.65 cm/ 
sec, a discharge temperature of 32. 3°C and a 16 meter constant depth 
region of interest. The total simulated time was 20. 6 hours. The 
results are summarized in Figures 18, 19 and 20. The surface 
velocities after 8.64 hours are shown in Figure 18, while Figure 19 
shows the velocities in a vertical plant (1=11) after 21.6 hours. 

Since the effects of wind were not included in this run and there 
were no effects of Jocassee, the velocities are by the discharge 
velocity and the buoyant plume. The isotherm comparison of the 
archival run and this simulation run is shown in Figure 20. The 
ambient temperature is 29°C. The archival isotherms are higher 
than the predicted isotherms. This is expected since this run 
(without wind) was not undertaken primarily for comparison pur- 
poses, and the simulation did not reach steady state. 

2. Run Number L002 

This run is similar to 'LOOT but with variable depth of the region 
of interest. The simulation time was only 8.64 hrs. This time is 
long enough to determine how well the model handles a variable 
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depth domain. The surface, velocities and vertical velocities (at 
1=11) are shown respectively in Figures 21 and 22. For complete-, 
ness, the simulated surface isotherms are compared with the archival 
isotherms. This comparison is shown in Figures 23. For the same 
reason described in the previous subsection, the comparison cannot 
be expected to be more precise. 

3. Run Number LOO 3 

Run L002 is now repeated with the effects of wind included. The 
total time simulated was only 8.00 hrs. The reason for the short 
simulation time is as explained above . The surface and 
vertical velocities (J =7) are shown in Figures 24 and 25 respectively. 
A comparison of Figures 21 and 24 shows the small effects of wind on 
the surface velocities after such a small simulation time. The wind 
speed is 4 m/sec. A comparison of the predicted and archival iso- 
therms is shown in Figure 26. Compared with the other isotherm 
comparisons, this figure appears to have improved the difference be- 
tween measured and predicted. 

Archival Runs and Results 


Two runs were made (LQ04 and L005) using the archival data base 
described in the last section. 

1. Boundary Conditions 

The following specifications were used at the main input boundaries. 

a. At 1=1, J=8 to 1 9 for all depths (K) 

The water velocity at the Jocassee end is specified as 6. 9 cm/ 
sec. Open boundary is specified for temperatures as this end. 

b. At 1=11, J=1, K=3 

This discharge velocity is specified as 5.65 cm /sec, and the 
dishcarge temperature also specified as 32.4°C. 

c. Finally, at 1=13, J=7, K=1 to 3 

Keowee dam discharge velocity is specified as 10.4 cm /sec. 

2. Run Number L004 

This is the first archival run and was made with variable depth 
topography for the near-field region of interest. 

3. Run Number L005 
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Run L004 is repeated but for constant depth region of interest. 

The domain was cut off at 16 meters from the water surface. This 
is the depth of the thermocline for Lake Keowee. 

4. Archival Results 

Figure 27 shows the comparison between the measured (archival 
> 9/10/75) and calculated (L004) isotherms* The ambient temperature 
was 29°C. This comparison can be seen to have improved when 
compared to those of the preliminary runs already described. The 
surface velocities are similar to 'L005 1 which are shown in Figure 28. 
Flow through the Keowee dam can be seen very clearly on the right 
boundary of the figure. Figure 29 shows the velocity field at J=7, 
which also shows the flow through the Keowee dam. The isotherm 
comparisons of measured (9/10/75) and predicted (L005) is shown 
in Figure 30. This is the best comparison of the runs described so 
far. The agreement is particularly good for the 29. 5 and 30°C iso- 
therms. The shapes of the measured and predicted 30. 5°C isotherms 
are similar. However, the measured area under the isotherm appears 
bigger. This discrepancy can be explained in terms of the nearness 
of this isotherm to the discharge. The total time simulated was 32.4 
hrs, which is also the time when steady state was reached. The 
temperature profiles for locations 1=11, J=7, K=1 to 4 through the 
simulated period are shown in Figure 31. The temperature profiles 
are vertical as expected since computations were carried out above 
the thermocline. A comparison of the surface velocities at this loca- 
tion is given in Figures 32. It can be seen from this figure that the 
u-component velocities predominate. It can also be seen from this 
figure that all the velocities stabilize after about 30 hours, showing 
steady state conditions. A similar plot. Figure 33, for the same 
variables but for location 1=11, J=2 and K=1 is shown. Similar con- 
clusions hold for this figure. 

L001 through L005 were archival runs. The following section dis- 
cusses the results of simulations for the August data base (L006) 
and February data base (L007). Both these runs were made in 
stabs of 24 hrs of simulation at a time for a total period of 48 hours 
of simulation in each case. Results were stored and printed at the 
end of every hour of real time. The results at the end of the first 
24 hrs simulation were used as initial conditions for the next 24 hrs. 

Summer Runs 

Simulation was started from 0000 hrs August 24. The model was 
started 'cold 1 i.e. , initial conditions used were zero velocities and con- 
stant temperature (29°C) throughout the domain. This conditions re- 
quired that the model had to be run for sometime before the effects of 
the cold start could be ignored. 

The inputs used in running the model have been discussed in the 
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previous section. The wind speed driving the upper layer varied from 
0.76 to 3.09 m/sec. The current or the flow through the Jocassee 
boundary varied from 0.61 to 1.12 cm /sec. The Oconee nuclear station 
discharge velocities and temperature were 7.42 cmVsec and 31.7°C re- 
spectively. (Average values were taken since the fluctuation was negli- 
gible.) A«i the above data were fed in every hour in running the model. 

The surface isotherms and velocities at each horizontal section 
(K=1 # 2, 3, 4 and 5) were plotted using a Calcomp plotter. Results at 
the end of each hour were stored and plotted. The results used for 
comparison were the surface isotherms plots. Measured surface isotherms 
were obtained from IR scanner digicolor maps. By drawing both the 
data and the predicted values on the same grid system direct comparison 
could be made. 

The first IR runs were made on August 24 from 0853 hrs EST to 
1002 hrs EST. The isotherms drawn for this data are shown in Figure 
34. The Isotherms shown are for 30. 5°C, 30. 0°C and 29. 5°C. Only the 
above three isotherms fell within the domain of interest and represent 
the water surface temperatures at 1002 hrs EST, August 24, 1978. The 
results produced by the model are shown in Figure 35 and 36. Figure 35 
shows the same isotherm i.e. , 30. 5°C, 30.0°C and 29. 5°C corresponding 
to 1007 hrs real time, it is seen that the temperatures predicted by the 
model are lower than the actual temperatures, showing that the plume 
spread was underpredicted. Figure 36 shows different predicted iso- 
therms corresponding to the same time. The temperatures of the iso- 
therms are 29. 5°C, 29. 3°C and 29. 1°C. Comparing Figure 36 with 34 
it is seen that 29. 5°C in the IR data agrees very well with the 29. 1°C 
isotherm predicted by the model. The measured 30. 5°C and 30.0°C iso- 
therms compared with the predicted 29. 5°C and 29. 3°C isotherms respec- 
tively. Hence, the errors in prediction in the three isotherms (30. 5°C, 

JiO. 0°C and 29. 5°C measured) are approximately 1°C, 0. 7°C and 0. 4°C 
respectively. 

The IR digicolor map for the August 24 afternoon data could not be 
used to construct isotherms because of bad resolution between different 
colors. Hence, the next useful comparison was made for the August 25 
morning data (0903-9553 hrs EST). The IR isotherms are shown in 
Figure 37. The same three temperatures, namely 30. 5°C, 30. 0°C and 
29. 5°C, are represented here since the other isotherms lie beyond the 
domain of interest. Figure 38 shows the same isotherms as predicted by 
the model. The spread of the 29. 5°C isotherm is overpredicted while 
the 30. 0°C isotherm is underpredicted. However, Figure 39 shows iso- 
therms for 29. 90°C, 29. 70°C and 29. 6°C as predicted by the model. 

These isotherms compare very well with the IR isotherms in Figure 37. 

The errors in this case are . 6°C for the 30. 5°C isotherm, 0. 3°C for 
the 30.0°C isotherm and -.1°C for the 29. 5°C isotherm. This shows a 
considerable improvement in the predictions as compared to August 24, 
showing the effects of cold start gradually vanishing. \ 
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Figure 40 shows the plot of the horizontal velocities of the surface 
as predicted by the model for 1007 hrs, August 24. During this time 
(1007 hrs EST), as can be seen from the Jocassee-pumped storage sta- 
tion flow and Keowee hydro station flow data r the flows through these 
stations were negligible. This leads to zero velocities at Jocassee boun- 
dary and zero velocity through the Keowee hydro discharge point in 
Figure 40. 

Figure 41 shows the surface horizontal velocities predicted by the 
model for August 25 (33.2 hrs run time). During this time (approxi- 
mately 1000 hrs, August 25) Jocassee-pumped storage station just started 
generating and Keowee hydro was generating. This accounts for the flows 
through the two boundaries. 

The main driving forces responsible for determining the shape of 
the isotherms are the ambient temperature, discharge temperature and 
flows through the Jocassee boundary. The wind is seen to affect the 
velocities in the upper layer only. This characteristic is displayed 
both in the data and the simulation results. Figure 34 (August 24 IR 
data) and Figure 37 (August 24 IR data) . Looking at the Jocassee and 
Keowee flows it is seen that the flow through Jocassee is negligible dur- 
ing both these periods. The flow through Keowee is negligible at 10.00 

а. m. on August 24 but is 9352 cfs on 10 a.m., August 25. The average 
discharge temperature is also lower on August 25 as compared to August 
24. This causes the area under the isotherms in Figure 37 to be lower 
than those in Figure 34. The same difference is seen between the cor- 
responding predicted isotherms, namely Figure 36 (August 24, 1978) and 
Figure 39 (August 25) . 

Winter Runs \ 

Simulation was started from (0000 hrs) February 24, 1979. The 
model was started using zero velocities and constant temperature (= 10°C) 
as initial conditions. 

The inputs used for this run (L007) have been discussed in detail 
in previous sections. The Oconee discharge velocity (average value) was 

б. 84 cm /sec and the discharge temperature (average value) was 18. 4°C. 
The wind speed varied from 1.61 to 4.52 m/sec and the Jocassee-pumped 
storage station boundary flow velocity ranged from 0.61 to 4. 14 cm /sec. 

Values of velocities and temperatures were printed at the end of 
every hour. Surface isotherm plots and velocity plots were generated 
from these results. Figure 42 shows the IR scanner isotherms for 1648- 
1651 hrs, February 27. The temperatures are 13. 0°C, 12. 5°C, 12.0°C, 
11.5°C and 11.0°C. The same isotherms as predicted by the model are 
shown in Figure 43. These correspond to 17.12 hrs after the cold 
start. Comparing these two figures it is seen that the isotherms pre- 
dicted by the model have a greater spread. Figure 44 shows predicted 
isotherms corresponding to 13. 75°C, 13.0°C, 12. 75°C, 12. 5°C and 12. 0°C. 


This figure agrees better with Figure 42 showing an error of prediction 
of 0. 75°C, 0. 5°C, 0. 75°C and 1°C respectively for each of the isotherms. 

Figure 45 shows IR scanner isotherms for February 28, 1979 corre- 
sponding to 0948-0957 hrs EST. The two isotherms within the domain 
are for 13. 0°C and 12. 5°C. Figure 46 shows the same two isotherms as 
predicted by the model. These correspond to 34.2 hrs of run time after 
the cold start at approximately 10 a.m., February 28. Figure 46 is in 
excellent agreement with Figure 45. Isotherms for IR data corresponding 
to February 28 afternoon could not be drawn because of bad resolution 
among different colors in the digicolor map. 

Figure 47 shows velocities at the surface as predicted by the model 
for February 27, 17.12 hrs after the cold start. The flows through the 
Jocassee and Keowee hydro boundaries agree with the data obtained. 
Figure 48 shows the same for February 28, 34.2 hrs after cold start. 
Excellent agreement is obtained in this case between the velocities shown 
on the map and measured flows at Keowee hydro station and Jocassee- 
pumped storage station. As in one of the summer runs, the isotherms 
are affected by the Oconee discharge, Keowee hydro and Jocassee-pumped 
storage station. Figures 42 and 45 (isotherms for February 27 and 28 
respectively) show that the isotherms have spread out further on the 
second date. Keowee hydro and Jocassee-pumped storage station flows 
are negligible on both the dates. The average Oconee discharge tempera^- 
ture during 1600 hrs, February 27, was lower than that during 1000 hrs, 
February 28. This accounts for the greater spread of the isotherms. 

The same is depicted in the results produced by the model. 
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Figure 1. Lake Keowee 






Figure 2. Coordinate and grid system 
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Figure 5. Lake Keowee (region of interest) showing inputs and outputs (for 3-D model) 










Figure 8. Flight plans for IR data (August and February missions) 
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Figure 10. Keowee hydro discharge data (August 25, 1978) 
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Figure 21. Velocities at K = 1, Lake Keowee (rigid-lid model) 
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Figure 22. Velocities at I - 11, Lake Keowee (rigid-lid model) 
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Figure 23. Isotherms: measured and predicted (temperatures above ambient) 
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RUN N0 2 LOO 3. 
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Figure 24. Velocities at K = 1, Lake Keowee (rigid-lid model) 
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Figure 25. Velocities at J =7, Lake Keowee (rigid-lid model) 
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Figure 26. Isotherms: measured and predicted (temperatures above ambient) 
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Figure 27. isotherms: measured and predicted (temperatures above ambient) 
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Figure 30. Isotherms: measured and predicted (temperatures above ambient) 


Reproduced from 
best available copy. 
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RUN HQ i L3Q 6. 

Q ISCHRRGE VELOCITY * 7.42CM/SEC 
DISCHARGE TEMPERATURE « 31.7®C 
HIND SPEED U1RX) * 3.Q9M/SEC 
CURRENT ( JQCRSEE FLOW 3 * 1 -t CM/SEC 
TQTRL SIMULATED TIME » I0.07HRS 



Figure 35. Isotherms at K = 1, Lake Keowee (rigid-lid model), 
simulations for August 24, 1978 
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RUN NO: L006 
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Figure 36. Isotherms at K = 1, Lake Keowee (rigid-lid model) 
simulations for August 24, 1978 
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Figure 40. Velocities at K =1, Lake Keowee (rigid-lid model), 
simulations for August 24, 1978 



RUN NG: LQG 6. 

DISCHARGE VELOCITY : 7.42CM/SEC 
OISCHRRGE TEMPERATURE t 31 . 7*0 
MIND SPEED CMRX) : 3.09M/SEC 

CURRENT! JQCASSE FLOW)* 1-1 CM/SEC 
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Figure 41. Velocities at K = 1, Lake Keowee (rigid-lid model), 
simulations for August 25, 1 978 







RUN NO: L007 
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Isotherms at K = 1, Lake Keowee (rigid-lid model) 
simulations for February 27, 1979 





RUN NO: L007 

DISCHARGE VELOCITY : 6. 84 CM /SEC 
DISCHARGE TEMPERATURE: 18.4°C 
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Isotherms at K = 1, Lake Keowee (rigid-lid model) 
simulations for February 27, 1979 








Figure 45. IR data corresponding to 0948-0957 hrs, February 28, 1979 
Lake Keowee 
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RUN NO: L 
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Figure 46. isotherms at K =1, Lake Keowee (rigid-lid model) 
simulations for February 28, 1979 
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Figure 47. Velocities at K = 1, Lake Keowee (rigid-lid model) 
simulations for February 27, 1979 
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Figure 48. Velocities at K = 1, Lake Keowee (rigid-lide model), 
simulations for February 28, 1979 
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Table 1. Monthly Gross Thermal Capacity Factors (Percent) 
for Oconee Nuclear Station 


MONTH 

1973 

1974 

1975 

1976 

1977 

Jan. 

0 

13 

37 

30 

90 

Feb. 

0 

28 

22 

49 

77 

Mar. 

0 

29 

57 

47 ■ 

72 

Apr. 

0 

26 

61 

27 

92 

May 

4 

5 

76 

33 

76 

Jun. 

9 

43 

78 

64 

46 

Jul. 

15 

38 

80 

59 

44 

Aug. 

15 

31 

82 

90 

36 

Sep. 

16 

46 

76 

77 


Oct. 

23 

r 

19 

75 

61 


Nov. 

23 

31 

92 

52 


Dec. 

27 

30 

92 

55 


Annual 

Average 

11 1 

28 

69 

59 



Data based on three-unit capability for entire period 

Gross Thermal A . . . nn 

Canacitv Factor = Actual MWH (t) x 100 

for P ONS(%) 7704 MW(t) x Hours in Period 












Table 2. Operating Conditions of Oconee Nuclear Power Plant 
at Lake Keowee (September 10, 1975) 


Flow Rate : 7087 m/min 

Inlet Temperature : 26.6 °C 

Outlet Temperature : 32. 4 °C 

The hourly averages of wind, air temperature and 
relative humidity: 


Time 

Wind 

Speed' (MPH) 

Wind 

Direction (°N) 


Relative 

Humidity 

9:00 

7.7 

15 

69.0 

100 

10:00 

6.8 

65 

74.0 

70 

11:00 

8.3 

80 

77.0 

60 

12:00 

7.2 

95 

77.0 

63 

13:00 

9.7 

115 

78.0 

66 

14:00 

11. 3 

205 

79.0 

64 

15:00 

- 

- 

77.0 

66 

16:00 

10.2 

75 

73.0 

84 








Table 3. Input Data for the 3-D Model (Lake Keowee) 


Flow Rate: 

118.12 m 3 /sec 

Inlet Temperature: 

25.6 °C 

Discharge Temperature: 

32.4 °C 

Discharge Velocity: 

5. 65 cm /sec 

Ambient Temperature: 

28.9 °G 

Depth at Discharge: 

30 meters 

Discharge Width: 

152. 4 meters 

Air Temperature: 

26.11 °C 

Wind Speed: 

12. 83 ft /sec 

Vertical Eddy Diffusivity: 

86.1 cm 2 /sec 

Horizontal Eddy Diffusivity: 

9. 2 x 10 3 /sec 



Table 4. Volume and Area Data for Lake Keowee 


WATER SURFACE 
ELEVATION 
(m) (ft.) 

SURFACE AREA 
(SQ. KM.) 

STORAGE VOLUME 
CIO 9 m 3 ) 

246.9 

— 

83.3 

1.42 

245.4 

■ 

78.8 

1.30 

243, 8 

800 

74.4 

1.18 

242.3 

■ 

70.1 

1.07 

240 .8 


65.8 

0.97 

239. 3 

785 

61.2 

0.87 

237.7 

780 

56.6 

0.78 

237.1 

778 

55.1 

0.75 

236.2 

775 

52.8 

0.70 

234.7 


48.9 

0.62 

231.6 


41.5 

0.48 

228.6 


35.3 

0.36 

225.6 

740 

27.8 

0.27 

222.5 

730 

22.5 

0.19 

219. 4 

720 

17.8 

0.13 

216.4 

710 

13.2 

0.08 

213.4 

700 

9.6 

0.05 

210.3 

690 

5.7 

0.02 

207.3 i 

680 

3.1 

0.01 

204.2 

670 

1.3 

0.004 

201.2 

660 

0.6 

0.001 

198.1 

650 

0.3 



* Full Pond 

** Maximum Allowable Drawdown 
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Table 5. Inflows ana OutfTowsr to Lake- Keowee, August 24, 1978 


TIME 

OCONEE 

DISCHARGE 

OCONEE 

DISCHARGE 

TEMP. 

NET JO CAS SEE 
FLOW 

KEOWEE 
HYDRO FLOW 

Aug. 24 
1978 

| (m 3 /min) 

°C 

(CFS) 

(CFS) 

12.00 AM 

8172. 9 

31.8 

-2638 

6636 

01.00 

8177.5 

31. 8 

-6059 

6696 

02.00 

8178.6 

31. 9 

-7982 

5416 

03.00 

8180.1 

32.1 

-12036 

5416 

04. 00 

8185. 4 

32. 3 

-12066 

6756 

05.00 

8179. 4 

32. 4 

-11985 

2276 

06.00 

8181.7 

32. 5 

-10598 

48 

07.00 

8313.0 

32. 3 

-5970 

48 

08,00 

8183.5 

32.4 

-995 

48 

09.00 

8176.4 

32.5 

100 

48 

10.00 

8183.9 

32.5 

100 

48 

11.00 

8177.1 

32.5 

100 

48 

12.00 PM 

8159.9 

32.6 

100 

48 

1.00 

8155. 9 

32.6 

100 

48 

2.00 

8247. 2 

32.7 

100 

48 

3.00 

8169.5 

32.7 

100 

48 

4.00 

8438.7 

32.4 

100 

48 

5.00 

7935.2 

30.7 

8726 

14224 

6.00 

7543.1 

31.0 

10905 

8432 

7.00 

8091.9 

30. 3 

9823 

4544 

8.00 

8119.6 

30. 2 

4557 

48 

9. 00 

8132. 8 

30. 2 

6482 

48 


8141. 9 

30. 2 

4183 

48 

11.00 

8143. 8 

30. 2 

100 

48 

12.00 AM 

8151.8 

30.3 

100 

48 














Table 6. Inflows and Outflows to Lake Keowee f ^August 25, 1978 


TIME 

OCONEE 

DISCHARGE 

OCONEE 

DISCHARGE 

TEMP. 

NET JOCASSEE 
FLOW 

KEOWEE 
HYDRO FLOW 

Aug. 25 
1978 

(m 3 /min) 

°C 

(CFS) 

f 

(CFS) 

12.00 AM 

8151.8 

30. 3 

100 

48 

1.00 

7952. 3 

28.9 

4938 

3576 

2.00 

7097.2 

218. 5 

838 

608 

3.00 

7098.4 

29.1 

100 

48 

4.00 

7098.8 

29.6 

100 

48 

5.00 

7103.7 

30.0 

100 

48 

6.00 

7102.6 

30.5 

100 

48 

7.00 

7104. 4 

30.7 

100 

48 

8.00 

7279.7 

30.9 

100 

3484 

9.00 

8012. 8 

30. 5 

100 

8624 

10.00 

7760.4 

30.7 

100 

9352 

11.00 

7759.2 

30. 9 

100 • 

1392 

12.00 PM 

8059.8 

30.7 

1 6 32 ,V 

2564 

1.00 

8130.6 

30.8 

5442 

8820 

2.00 

8112. 8 

30.9 

11553 

9500 

3.00 

8141.9 

31.1 

19048 

7624 

4.00 

8102. 2 

’ 31.3 

17180 

7856 

5.00 

8108.6 

31. 3 

10615 

8852 

6.00 

8107.5 

31.4 

10158 

8172 

7.00 

8110.1 

31.5 

10076 

2660 

8.00 

8114.7 

31.5 

21 85 

48 

9.00 

8126. 8 

31.5 

29. 83 

48 

10.00 

8148. 3 

31.5 

3151 

48 

11.00 

8141.9 

30.4 

2647 

48 

12.00 AM 

8141. 9 

30. 3 

100 

48 


Table 7. Meteorological Data for Lake Keowee, August 24, 1978 


Time 

(Hrs. from 
midnieht) 

Wind Speed 
(in/s) 

Air Temp 
C°C) 

Dew Pt. 
Temp. 
(°C) 

Solar 

Radiation 

(w/m 2 ) 

Wind 

Direction 

1 

1.788 

21. 9 

21. 94 

0.0 

30.0 

2 

2.76 

21.4 

21.39 

0.0 

10 

3 

1. 96 

20. 8 

20.83 

0.0 

30.0 

4 

2.53 

20. 8 

20. 83 

0.0 

20 

5 

7.45 

20. 3 

20.28 

0.0 

35.0 

6 

2. 27 

20.0 

20.0 

0.0 

20.0 

7 

1.38 

IS. 4 

19.4 

34. 85 

30.0 

8 

2. 85 

20.4 

20.0 

170.77 

15.0 

9 

1.56 

21.4 

i 

20.83 

326.43 

10.0 

10 

1.29 

23.3 

21.67 

592.45 

40.0 

11 

1.69 

26.1 

23.06 

710. 94 

20.0 

12 

i 

1. 96 

27. 2 

23.89 

815. 49 

30.0 

13 

2. 90 

28.1 

24. 44 

864.28 

40.0 

14 

2. 86 

28. 9 

25.56 

836.4 

15.0 

15 

2.90 

29.4 

25. 83 

752.76 

85.0 

16 

1. 82 

29.4 

25. 28 

599.42 

60.0 

17 

1.78 

29. 4 

25. 28 

397. 29 

70.0 

18 

1.42 

29.4 

26.67 

214. 9 

85.0 

19 

1. 47 

29. 9 

26.53 

81. 32 

75.0 

20 

‘i 

27.2 

25. 97 

5.81 

90.0 

21 

0.71 

25. 8 

25.28 

0.0 

30.0 

22 

1.82 

25. 28 

25.0 

0.0 

20.0 

23 

1.16 

24.17 

24.17 

0.0 

5.0 

24 

1.61 

23.89 

23.89 

0.0 

25.0 
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Table 8. Meteorological Data for Lake Keowee, August 25, 1978 


Wind Speed 
(m / s) 

Air Temp 
(°C) 

1.028 

23.1 

2.951 

23.1 

2. 414 

22.? 

2. 995 

22 2 

1.833 

21.9 

1. 833 

21.7 

2. 235 

21.4 

1. 922 

21.7 

1. 833 

23. 3 

0.760 

25. 3 

1.162 

26.9 

1.565 

28.6 

1.654 

29.7 

1. 967 

30.6 

2. 503 

30.8 

2.012 

30.8 

2.056 

30.8 

1.743 

30.0 

1. 967 

29.4 

1.654 

28.3 

2. 414 

27.2 

1. 341 

25.8 

1.073 

25.6 

3.085 

25.6 



Solar 
Radiation 
(w /m 2 ) 


23.0 

23.0 
22. 2 
22 . 2 

21.9 

21.7 
21.4 

21.7 

20.0 
17.2 
18. 3 

16.1 

16.9 
16. 7 

16.1 

16.1 

16.7 

20.0 

22.7 

21.1 

25.6 

25.6 

25.6 

25.6 









Table 9. Inflows and Outflows to take Keowee, February 27, 1979 



12.00 PM 


3 

4 

5 

6 

7 

8 
9 
1 
1 

12. 00 AM 


OCONEE 

DISCHARGE 

(nr' /min) 

OCONEE 

DISCHARGE 

TEMP. 

C°C) 

NET JOCASSEE 
FLOW 

(C.F.S.) 

KEOWEE HYDRO 
FLOW 

(C.F.S.) 

7505. 3 

18.6 

-14395 

48 

7498.1* 

18.5 

-18754 

48 

7492.0 

18. 4 

-18805 

48 

7492.0 

1 8. 5 

-18713 

48 j 

7491.6 

18. 3 

-18698 

48 

7494. 3 

18. 3 

-18688 

48 

7488. 2 

18.3 

-15939 

48 

7481. 8 

18.2 

3484 

3668 

7485.6 

18. 3 

16823 

17540 

7488.2 

18. 2 

13503 

8488 

7497. 7 

18. 3 

5470 

8096 

7504.1 

18. 3 

100 

2680 

7503. 4 

18.4 

100 

48 

7506.0 

18.5 

100 

i 

48 

7506. 4 

18.5 

100 

48 

7503. 4 

18.5 

100 

48 

7501. 9 

18. 4 

100 

48 

7507.5 

' 18.4 

100 

48 

7511.0 

18.4 

100 

48 

7516. 2 

18.4 

100 

48 

7518. 9 

18.3 

100 

48 

7520. 4 

18.3 

100 

48 

7516.6 

18.2 

100 

48 

7509.4 

18.2 

100 

48 

7507. 2 

18.2 

-4382 

48 










Table 10. Inflows and Outflows to Lake Keowee, February 28, 1979 


OCONEE 

DISCHARGE 

OCONEE 

DISCHARGE 

TEMP. 

NET JOCASSEE 
FLOW 

KEOWEE HYDRO 
FLOW 

(m 3 /rain) 

(°C) 

(C.F.S.) 

(C.F.S.) 

7507. 2 

18. 2 

-4382 

48 

7616. 4 

19.8 

-18050 

48 

7655.5 

19.5 

-18621 

48 

7356. 3 

19. 8 

-18602 

48 

7233.4 

20.1 

-18692 

48 

7494.7 

19.5 

-18559 

48 

7605.3 

19.1 

-15139 

48 

7618.6 

19.1 

3045 

48 

7300. 2 

19.5 

19245 

5788 

7239.8 

19.5 

8338 

48 

7241. 9 

19. 4 

100 

48 

7243.0 

19. 3 

100 

48 

7246. 2 

19.1 

100 

48 

7239.0 

19.0 

100 

48 

7241.0 

18.9 

100 

48 

7235. 3 

18. 9 

100 

48 

7207. 8 

18. 9 

100 

48 

7196. 3 

18.7 

100 

48 

7200.0 

18.9 

100 

48 

7211.5 

19. 2 

100 

48 

7212. 

19.0 

100 

48 

6829. 4 

18. 3 

100 

48 

7189.1 

18.5 

100 

48 

7287. 8 

18.1 

-1544 

48 

1 70 36.7 

• SJ ■ 

15. 8 

-15141 

48 
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Table 11. Meteorological Data for Lake Keowee, February 27, 1979 


Time 

(Hrs from 
midnight) 

Wind Speed 
(m /s) 

Air Temp 
(°G) 

Dew Pt. 
Temp. 
t°C) 

Solar 
Radiation 
(w /m 2 ) 

Wind 

Direction 

(Degrees) 

1 

1.833 

- 0.33 

- 2.78 

0.0 

15 ° 

2 

1.073 

- 0.72 

- 1.67 

0.0 

75 ° 

3 

2 . 325 

- 1.61 

- 1.61 

0.0 

60 ° 

4 

. 1.565 

- 2.22 

- 2.28 

0.0 

15 ° 

5 

2.056 

- 1.83 

- 1 . 89 

0.0 

50 ° 

6 

1.788 

- 2.17 

- 2.22 

0.0 

85 ° 

7 

2.012 

- 2.72 

- 2.78 

20 . 94 

85 ° 

8 

2 . 280 

- 1.67 

- 2.78 

195 . 39 

60 ° 

9 

0.626 

0.01 

- 3 . 33 

369 . 85 

5 ° 

10 

1 . 386 

3.06 

- 2.22 

544. 31 

75 ° 

11 

1.609 

5.83 

- 2.22 

655 . 96 

15 ° 

12 

1.788 

8.83 

- 1 . 39 

725.75 

40 ° 

13 

3,129 

11.06 

- 2.78 

746 . 68 

80 ° 

14 

2.593 

12.28 

- 5.0 

704 . 81 

70 ° 

15 

1.520 

13.39 

- 5.56 

579 . 20 

80 ° 

16 

1.207 

13.89 

- 5.56 

383 . 81 

75 ° 

17 

1.565 

13 . 83 

- 5.61 

146.55 

55 ° 

18 

1.609 

13.72 

- 3 . 33 

20 . 94 

15 ° 

19 

2.056 

11.72 

- 4.44 

0.0 

30 ° 

20 

1.162 

r ,72 

- 2.78 

0.0 

25 ° 

21 

2 . 772 

0 , 33 

5.28 

0.0 

55° 

22 

2 . 861 

7.78 

5.56 

0.0 

55 ° 

23 

2 . 995 

7.00 

5.28 

0.0 

50 ° 

24 

1 . 386 

5.28 

3 . 89 

0.0 

60 ° 




r Lake Keowee, February 28, 1979 


Dew Pt 


4.17 
4.17 

2.22 
1.67 
2 . 22 
3.06 
0.28 
-2.78 
-3. 89 
-5.00 
-5.00 
-3.89 
-3. 89 
-3. 89 
-3.33 
-2. 78 
-1.67 
-0.56 
1. 39 
2.28 

2.22 

2.22 

2.22 

2.22 


Solar 
Radiation 
(w An 2 ) 

Wind 

Direction 

(Degrees) 

0.0 

55 

0.0 

55 

0.0 

60 

0.0 

50 

0.0 

80 

0.0 

75 

20. 94 

60 

153.52 

60 

279.13 

85 

439.64 

60 

676. 90 

70 

907.18 

45 

537. 33 

60 

565.24 

50 

690. 85 

45 

0.0 

45 

0.0 

45 

0.0 

5 

0.0 

65 

0.0 

55 

0.0 

85 

0.0 

50 

0.0 

65 

0.0 

70 
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